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ABSTRACT 

In light of recent indications for a large population of obscured active galactic nuclei (AGNs) , we revisit 
the mean radiation efficiency from accretion onto supermassive black holes (SMBHs), (e), applying a 
bayesian approach. We use the integrated comoving energy density emitted by AGNs and compare 
it to the mass density of local SMBHs. When considering only optically-selected unobscured AGNs, 
we derive logio[(e)]=-1.77lo or (e) = 1.7to;4%. Using the AGNs selected using hard X-rays, which 
include unabsorbed and Compton-thin absorbed AGNs, we find logio[(e)]=-1.20lo'i[Q or (e}=6.4l^ 3%. 
Using optically-selected AGNs, and correcting for the obscured population, we inferr logio[(e)]=:- 
1.171 

0.08 (^)— which we consider our best estimate. We also repeat this calculation for 
intrinsically luminous AGNs (Mb < —23, quasars), comparing to the SMBH mass density in local 
elliptical galaxies, and find logio[(e)]=:-1.27lQ or (e)=5.4j;^'2%. We discuss how our results can be 
used to estimate the mean spin of accreting SMBHs, both assuming the standard thin-disk model of 
accretion onto black holes and also comparing to more recent simulations. Our results can rule out 
maximally rotating SMBHs ((a) = 0.998 '^^' ) at the >98% confidence level, as well as high rotation 
values ((a) > 0.9^?*-) with >92% confidence. Our preferred values of (a) are ~0. 25-0.60^?*, ahhough 
even these might be overestimated. Hence, we find that on average, SMBHs are not rapidly spinning 
during accretion. Finally, using an independent measurement of Eddington ratios, we estimate the 
SMBH e-folding time for the brightest AGNs (quasars) to be (r) = lOOtgQ^ Myr. 

Subject headings: galaxies : active -galaxies: nuclei ~ quasars: general -black hole physics - cosmology: 
miscellaneous - 



1. INTRODUCTION 

Active galactic nuclei (AGNs) are believed to be 
powered by accretion onto supermassive black holes 
(SMBHs) , and are expected to leave beh i nd dorman t 
SMBHs (|Salpete j IT961 iLvnden-Bell IT961 iReei fl98l . 
The observed energy per comoving volume radiated by 
AGNs, the density of local SMBHs and the radiation 
efficiency of accreting SMBHs can all be related, using 
co nservation of e nergy, in an elegant argument proposed 
bv lSoltanI ([T98l ). 

The original form of the "Soltan argument" esti- 
mated the density of relic SMBHs from the distribu- 
tion of observed quasar fluxes (magnitudes) and mak- 
ing an educated guess for the efficiency. Nowadays, 
SMBHs are believed to reside at the centres of galax- 
ies, and estimates of the SMBHs masses are possi- 
ble ( e.g. Kormendv fc Rich stond [19951 lEckart fc Genzell 
119961 : 'Magorrian ct al. 199^. More accurate measure- 
ments of the AGN luminosity functions (LFs) exist, 
meaning that the least constrained parameter is the ra- 
diation efficiency. Using variations of the Soltan argu- 
ment, the efficiency can be estimated from the com- 
parison between the estimated total radiated energy 
from AGNs with the local SMBH (e .g. Elvis ct al 2002; 
Yu fc Tremain3l2002l : iMarconi et"al][2004 iShankar et all 
20041) . 

The discovery of a significant population of obscured 



AGNs at a range of redshifts (see Section |4|) has impor- 
tant consequences for the global energetics of quasars and 
accretion onto SMBHs, given that the total energy out- 
put from the obscured population is comparable to that 
of the unobscured population (if not greater) . 

In this article, we address the issue of how this pop- 
ulation can have a significant effect on the mean radia- 
tion efficiency from accretion onto SMBHs. We apply a 
bayesian inference technique to derive posterior probabil- 
ity distribution functions (PDFs) for the mean radiation 
efficiency of SMBHs, (e). 

In Section [2] we outline how the total radiated energy 
density and the local SMBHs density can be used to es- 
timate (e). In Section [3] we discuss the relative contribu- 
tions of AGNs at different redshifts to the total radiated 
energy density. Section |4] summarises the current knowl- 
edge on the obscured AGN population. In Section [5] we 
explain in detail how we inferr (e), using a bayesian ap- 
proach, and how we combine different estimates of (e). 
Section [S] discusses how the mean SMBH spin during pe- 
riods of accretion can be estimated from (e) . Finally, in 
Section[7]we estimate the e-folding time for quasars, and 
in Section [5] we summarise our results. 

Throughout this article we refer to powerful AGNs as 
quasars where, for unobscured AGNs, the dividing line is 
generally taken at Mb = -23 [L^^ =7 x lO^^ Hz^i). 
We assume a ACDM cosmology with the following pa- 
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rameters: h = iJo/(100 km s"^ Mpc^^) = 0.7; = 0.3; 
^ 0.7. 

2. THE EFFICIENCY OF ACCRETION ONTO SMBHS 

We outline here the main steps of the Soltan argument, 
in the form that we will use to inferr the mean radiation 
efficiency onto SMBHs, (e). 

The total (comoving) energy density radiated by AGNs 
can be related to the integrated local SMBH density, 
p»,iocj using three assumptions: i) that SMBHs grew all 
their mass by accretion, thus neglecting mergers between 
black holes (or rather taking the simplifying assumption 
that mergers conserve SMBH mass) ii) that the mean 
efficiency (e) is constant with luminosity and over the 
Hubble time and iii) that AGNs are the manifestation of 
these periods of growth by accretion. 

For an individual SMBH, accretion with radiation ef- 
ficiency e will lead to a bolometric energy: 

/dm T , 
e— c^dt, (1) 

while a relic SMBH will be left behind, with mass, 

^ dm , 

{l-e)—dt. (2) 

The radiated energy no longer contributes to the accreted 
mass energy, which the distant observer will see as the 
SMBH mass, hence the (1— e) term. Rather than tracking 
individual sources, we are interested in the relic SMBH 
mass density (p,,acc) inferred from the summed comov- 
ing energy density emitted by AGNs over the age of the 
Universe, which can be estimated from the AGN LF. 
Converting to observable quantities, we obtain: 

P.,acc = ^^^^y§^ I J C^vL, <f>{L,, z)dL,, (3) 

where z) is the LF for all AGNs, Ci, is the bolomet- 

ric correction and (e) is the mean efficiency. The bolo- 
metric correction (defined so Lboi — CuvL^) is necessary 
to convert from the monochromatic LF to the bolometric 
luminosity density. 

If only an unobscured AGN LF is available, one can 
instead correct for the missing obscured AGNs: 

P..acc = -^^^^y^(l+(i?)) j ^dz j CliyL,ML.,z)dL, 

(4) 

Here (pn{L^,z) is the unobscured AGN LF. (R) is the 
mean ratio of obscured to unobscured AGNs so the term 
{1+{R)) is a parametrised way of accounting for the ob- 
scured population. Estimates exist for the parameters 
Ci, and (R). As mentioned earlier, the least constrained 
parameter is (e), so we treat it here as a free parameter. 

3. RADIATED ENERGY DENSITY 

A crucial quantity for the estimation of (e) is the 
total energy density radiated by AGNs, U (the time- 
integrated luminosity density). The AGN population 
shows a strong evolution so that the contribution of the 
low-redshift AGNs {z ^0.3) to the energy density is very 
modest. 

Amongst optically-selected AGNs, the luminosity den- 
sity of AGNs strongly incre ases with red shift, with a peak 
around z ^ 2-2.5 (e.g. Wolf et al.ll2003D . and for the en- 
tire AGN population over cosmic time, U is dominated 



by the AGNs at z =1-3 (particularly those around L*) 
with the lower-redshift AGNs contributing a smaller frac- 
tion. 

This is illustrated in the left panel of Figure [1] which 
shows the cumulative energy density U (z) [relative to 
U{z=0)] for several optical LFs and one hard X-ray LF. 
For the optical surveys, all defined in the B-band (centred 
on 4400 A, or 2.8 eV), U{z) is defined as: 

/""'"^ dt f 

Ub{z) = / dz''^^' / '-^Bt^B^^B (f>B{Li^B,z')dL^B- 

(5) 

where hpi^B —2.8 eV {hp is Planck's constant). L^b, 0b 
and Cb are, respectively, the luminosity density, LF and 
bolometric correction at the B-band. Plotting U{z)/U{0) 
removes the dependence on the bolometric correction and 
shows the relative contribution at different epochs. 

As can be seen in the left panel of Figure [H for 
optically-selected AGNs, approximately 70% of U{z=0) 
is radiated between z = 1 and 3, with 20% radiated at 
z < 1 and 10% at z > 3. 

However, similar studies undertaken at hard X-ray en- 
ergies (in the 2-10 keV band) have shown that the lower- 
luminosity AGNs evolve differently to high-luminosity 
AGNs, with their activity peaking at lower redshifts (e.g. 
iUeda et al.ll2003l). A gain, we illustrate the effect by plot- 
ting U{z) (Figure [U right panel), with U{z) defined as: 

f^'^"" dt r 

Uxiz)^ J ^ d^dz'y CxLxMLx,z')dLx. (6) 

where Lx, 4>x and Cx are the hard X-ray luminosity, 
LF and bolometric correction, but note that here iboi — 
CxLx- 

From the right panel of Figure [T] we can see that only 
about 50% of U was radiated between z = 1 and 3, with 
~45% radiated at z < 1 and ~5% at z > 3. 

Thus, from optically-selected LFs one expects the 
AGNs at 1 < z < 3 to dominate, although the hard 
X-ray LF indicates that the contribution of z < 1 AGNs 
is similarly important and should not be overlooked. 

4. DEMOGRAPHY OF OBSCURED AGNS 

Obscured AGNs are believed to be accreting SMBHs 
where, between X-ray and near-infrared wavelengths, gas 
and dust block the line of sight to the central region. The 
strong continuum from the accretion disk and the broad 
emission lines suffer from heavy extinction by dust, while 
the narrow emission lines may sometimes be visible, de- 
pending on the obscuring geometry. Photoelectric ab- 
sorption by the intervening gas and dust will supress the 
soft X-ray emission (^5 keV). If the column density is 
large enough to be optically-thick to Compton scatter- 
ing^, then even the hard X-rays (^5 keV) will be heavily 
affected. 

The obscured AGNs are indistinguishable from galax- 
ies in optical and near-infrared imaging surveys, which 
has made their identification significantly harder than 
that of their unobscured counterparts. There are mainly 

^ A source is Compton-thick if the absorbing column density, 
^Hi is > 1/t =1.5 X 10^* nr"'^, where aj' is the Thomson electron 
scattering cross-section. X-ray absorbed, Compton-thin sources 
have 10^ ^ Nn ^ 10^* m~^, and X-ray unabsorbed sources 

have A^H ^ lO^^ m^^ 
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four methods for fooking for these objects, which we sum- 
marise in chronofogical order, and we concentrate here on 
the most recent results. 

Radio selection (wavelength ^6 cm, energy ^20 /zeV). 
Dust and gas are transparent to moderately high- 
frequency radio waves, so radio-loud obscured AGNs 
can be readily identified, in the form of narrow-line 
radio galaxies. Amongst the radio-loud population, 
the obscured to un obscured ratio is found to be ^1:1 
(|Willott et al.l l2000f ). However, the fraction of unob- 
scured AGNs tha t are radio loud is small (~10%, e.g. 
Ilvezic et al.ll2002[ ). so that the narrow-line radio galaxies 
are expected to represent a similarly smal l fraction of t he 
obscured AGN population fe.g. [Zakamsk a~e"t al.ll200't ). 

Selection at hard X-ray energies (^2-10 keV). This 
method is sensitive to the Compton-thin absorbed AGN 
population, and it has iden tified large numbers of sou rces 
in deep X- ray surveys fseelBrandt fc Hasingeij[2005l . for 
a review). iGilli et all (|2007l ) find the ratio of Compton- 
thin absorbed to imabsorbed AGNs to vary from ~1:1 at 
the high luminosity end (Lx~ lO^'^ Lq) to ~4:1 at the 
low luminosity end (Lx~ 10^ L©). Using the hard X-ray 
background (at energies '^3-100 keV) as an extra con- 
straint, these authors find that the Compton-thick and 
Compton-thin absorbed populations are likely to have 
similar numbers (i.e ^--^2:1 at high luminosities and up to 
~8:1 at low luminosities). 

Selection using optical spectroscopy (wavelengths 
~3500-7500 A, energies ^1.5-3.5 eV). If the geometry 
of the dust is such that the narrow emission lines are not 
obscured, these are detectable via optical spectroscopy. 
Using the spectr oscopic database from the Sloan Dig- 
ital Sky Survey, IZakamska et al.l (|2003f ) were able to 
select large numbers of obscured luminous AGNs (ob- 
scured quasars). Using the [O III] 5007 A line strength 
as a isotropic proxy for the intrinsic AGN luminosity, 
the ratio of spectroscopically-selected obscured quasars 
to unobscured quasars is ~1. 2-1.5:1 (|Reves et alJl2008l 
for L[o HI] ^ 10^ Lq). At lower [O III] luminosities 
(-^[O III] ~ 10^ L©) ISimpsonI (j2005[ ) estimate the ob- 
scured to unobscured ratio to be ~4:1. 

Mid-infrared selection (wavelengths ^3-30 /im, en- 
ergies 0.04-0.4 eV). The extinction due to dust be- 
comes very small at these wavelengths, so mid-infrared 
selection should be sensitive to the radio-quiet popu- 
lation, including obscured AGNs showing no narrow 
lines. Recent work has suggested i t is also sensitive to 
the Compton-thick population (e.g. iPoUetta et al.|[2006t 
[Martfncz-Sansigre ct al. 20Q3)- However, care must 
be taken to separate starburst contaminants from real 
AGNs, since these star-forming galaxies also show bright 
mid-infrare d emission. Combining mid-inf rared and ra- 
dio criteria, iMartmez-Sansigre et al] (|2008D found a -2:1 
ratio at the high-luminosi ty end {fL,^ s ^tm 10^^ L©), 
while [Ponlev et al] ()2005f ) found a ~4:1 ratio amongst 
lower luminosity AGNs (i^L^ 24 nm ~ 10^° L©). 

B ased on mid-infra red selection with no radio crite- 
ria, |Lacx£tIS] l|2007| ) found a —2:1 ratio amongst bright 
AGNs (with a median valu e i/L^, 5 um ~ 5 x 10^^ L©) 
Using a power-law criterion. I Alonso-Herrero et al.l ()2006f ) 
also found a ~2:1 ratio amongst powerful AGNs (70% of 
their s ources having Lir > 10^^ L©) while iDonlev et all 
()2007f ) find a 4:1 ratio for less luminous sources (with 



typically Lx ^ 5 x 10^ L© or expected total infrared 
luminosities Lir ^ 10^^ L©) but only 2:1 if the sources 
are required to be X-ray detected. 

A possible picture fitting these observations is the 
following: amongst luminous AGNs (quasars), the ob- 
scured to unobscured ratio seems to be —2:1, with 
about half of the obs c ured quasars being Co r npton - 
thick (|Gim et all 120071 : iMartfnez-Sansigre et al] 120071 ). 
and about half showing no narrow emission lines 
(|Martfnez-Sansigre et al.ll2008l which explains the —1:1 
ratio found from spectroscopically-selected samples by 
Reyes et al., 2008). Amongst radio-loud quasars, the 
obscured to unobscured ratio is —1:1 with almost all 
obscured quasars (narrow-line radio galaxies) showing 
narrow emission lines. An explanation for this is that 
the lack of narrow lines is due to kpc-scale d ust (e.g. 
IMartfnez-Sansigre et al]|2006l: iRigbv et al.ll2006D and ra- 



dio jets are efficien t at clearing dust (|Baker et al.l 120021 : 
IWillott et al.ll200^ . so the fraction of radio galaxies ob- 
scured by kpc-scale dust is small. 

Amongst the low-luminosity AGNs, from X-ray, opti- 
cal spectroscopy and mid-infrared studies the obscured 
to unobscured ratio seems to be —4:1, c onsistent with 
the r atio inferred in the local Universe (e.g. lRisaliti et al"] 
[1999"). However, modelling of the hard X-ray background 
allows this ratio to be as high as 8:1 provide d that half 
of th e absorbed AGNs are Compton-thick (|Gilli et al.l 
[2007b . 

All the evidence points towards an obscured frac- 
tion that increases with decreasing AGN lu minosity 
(whic h was noted a long time ago, see iLawrence fc Elvig 
Il982f ). Whether the obscured fracti on evolves or not 
with redshift is still unclear (see e.g. lUeda et all [20031 : 
iTreister fc U77^ [200fil : iGilh et all [20f)l . The studies 
cited above range from z < 0.3 to 2; — 2j^and luminous 
AGNs are present in t h e entire range ([Za kamska ct aO 
120031 iLacv et al.l l2007l : IMartmez-Sansigre et all l200i l. 
Lower-luminosity AGNs, however, are typically only 
found at low and intermediate redshifts {z ^1), as ex- 
pected from fiux-limited samples. Thus, it is difficult at 
this stage to judge whether the obscured fraction varies 
with redshift. 

5. INFERRING (e) 

Having discussed the radiated energy density from 
AGNs and the inferred ratio of obscured to unobscured 
sources, we proceed to use these to estimate the mean ra- 
diation efficiency of AGNs. We will inferr different values 
of (e), from the optical LF, with and without adding ob- 
scured AGNs, as well as the hard X-ray LF. This will 
illustrate the effect on (e) of the obscured AGN popula- 
tion. 

Under the assumption that the local SMBHs grew by 
accretion: 



P.,loc 



(7) 



With p.^acc given by Equations [3| or [H we can esti- 
mate the value of (e) which yields the closest agree- 
ment between p.joc and p,,acc- For convenience we define 
u = U/C,^, where U is defined by Equation [5| or [6l This 
is done to keep separate and explicit the uncertainties in 
the LF (represented by u) and on the bolometric correc- 
tion {C„). 
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We now apply Bayes' theorem, to obtain the posterior 
PDF of (e), given the measurement of /3»,ioc in the local 
Universe and marginalising over the parameters u, (R) 
and Ci,: 



gaussian priors around these values, with the former hav- 



ing typical uncertainties 



=0.05 and 0.10 for the 



B-band and hard X-rays, respectively, and the latter 
=0.08 and 0.11 for the same bands. 



P((f)|p..ioc) = f^^P(p.,ioc| (e)) 



P(P.,loc) 

P((^)) 



P(u)du X 



P(p.,loc) 

P((i?))d(i?) / P(a)P(p.,ioc| {e),u,{R),C,)dC,. (8) 



For our purpose, p.joc is our data, so that P(p«,ioc) 
is the evidence, P(m), P{{R)) and P{C,^) are the prior 
PDFs for u, (R) and Ci, respectively. We have made 
use of the fact that our parameters are independent, 
so that P((i?)| (e))=P((i?)), P(u| (e))=P(u) and P(C^| 
(e))=p(C^). P(p.joc| (e), w, (i?), C^) is the likelihood of 
P.joc given (e), {R) and C^: 



1 ( P».loc~P«,a 



P(p.,ioc| (e), M, (R), G^) = -77r-W2 

From Equation |4] and the definition of u, p,,acc is: 



(9) 



p. 



(1 - (e)) 



(6) C2 



(1 + 



(10) 



(11) 



depending on whether obscured AGNs are to be added 
on via the {1+{R)) term. Note that we have made the 
simplifying assumption that (R) is independent of lumi- 
nosity, this will be addressed in Section [STU 

5.1. Priors 

We now assign a prior PDF to each of the parameters 
Ci/, (R), u, and (e), considering three scenarios: a) con- 
sidering only unobscured AGNs b) considering the X-ray 
absorbed and unabsorbed AGNs, but not the Compton- 
thick ones and c) attempting to include all obscured and 
unobscured AGNs d) same as c) but considering only the 
bright AGNs (quasars). 

5.1.1. P(a) 

In this work we are interested in the bolometric output 
from accretion onto the SMBH, which we assume pro- 
duces X-ray, optical and ultraviolet light predominantly 
and emits isotropically. We also assume that the dust 
responsible for the observed infrared emission (and the 
orientation-dependent SED) is on a larger scale, and not 
directly relevant to the accretion disk. Thus, the bolo- 
metric corrections for accretion onto the SMBH must 
convert from B-band luminosity to bolometric luminosity 
integrated over the optical, ultraviolet and X-ray range 
only 

We ado pt the lumino sity-depen dent bolometric correc- 
tions of iMarconi et al.l |2004. ) and lHopkins et all (|2007D . 
both of which follow these assumptions. We assume 



5.1.2. P(m) 

We estimate the total radiated energy density from 
unobscured AGNs by using the B-band LF, and the total 
radiated energy density from unabsorbed and absorbed 
(but Compton-thin) AGNs using the hard X-ray LF. The 
evolution of the LF is such that the luminosity density at 
z > 5 is similar to that at z < 1, but the lookback time 
is much shorter, meaning the energy emitted at z > 5 is 
negligible. In other words, the integral in Equation |4] has 
converged by z = 5. However, the LFs are not as well 
constrained in the range 3 < z < 5, and the functional 
form chosen to represent the evolution (e.g. a 2nd or 
3rd order polynomial) can have unphysical effects. We 
therefore prefer to restrict the range over which the LFs 
are integrated to < z < 3. The missing radiated energy 
is expected to be ^10% (see Section and therefore 
negligible compared to our uncertainties. 

Due to the fact that the bolometric corrections are 
luminosity-dependent and need to be included inside the 
integral when calculating U , we quote here the product 
of u and assuming the mean value of . The uncer- 
tainty quoted is only due to u, estimated by integrating 
different LFs, while the uncertainty in Ci, is encoded in 
the prior P{C^). 

For the B-band, due to the differe nt surve y areas and 
sensitivities, we use the LFs of Bovle et al.l (l200ft) and 
iCroom et al.l (pOOl ) at z < 1 and of lWolf et al.l (|2003h 
in the range 1 < z < 3 (using both pure density and 
pure luminosity evolution). The LF is integrated for 
aU magnitudes brighter than Mb < -18 (L^b>8x10^° 
W Hz^^). The estimated total energy from the B-band 
LFs is found to be MC„=(l.lib0. 2)xlO'^^ J Mpcy^ us - 
ing the bolometric corrections of IMarconi et al.l (|2004l ). 
and uC^=(1.6±0.4)x lO^^ J Mpc'^ when using from 
iHopkins et all (|2007D . 

If only the most powerful AGNs are considered 
(quasars, brighter than A/b=-23 or LyB=7xlO^^ 
W Hz -^, we find mC'^=(0 .9±0.2) x 10^^ J Mpc'^ (C^ 
from IMarconi et ail 120041) and MC^= (1.4ib0.3)xlO'^^ 
J Mpc-3 (C^ from 'Hopkins eLaljll^. 

We note all of the B-band LFs are likely to be missing 
the population of reddened AGNs. These sources suffer 
from moderate amounts of extinction in the visual band. 
Ay ^1 — 5, which suppresses the AGN signatures at 
optical but not at near- infrared wavelengths ('^1-2 /im, 
energies ~0.5-l eV). Brown et al. (2006) have estimated 
these to comprise ~20% of the unobscured population, 
so the above estimates of the energy density include this 
uncertainty. 

For the h ard X -ray L F', we use the L Fs o f 
lUeda et al.1 (|2003D and iSilverman eTal! (|2008D . 
The LF is then integrated over Lx >3 x lO^'' W. 
The estimated radiated energy density, which in- 
cludes Compton-thin absorbed AGNs, is then 
uC„=(4.5ib0.1)xlO ^^ J Mpc-3 (using from 
IMarconi et 311120041) and MC,.=f6.4±0.1) x 10^^ J Mpc^^ 
(C„ from IHopkins et al.ll2007h . 
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5.1.3. P((i?)) 

As we have discussed in Section 21 the value of (i?) 
inferred amongst high-luminosity AGNs is ~2, while 
amongst low-luminosity AGNs it is ~4. In Equations |4] 
and [Til we have made the simplification of keeping (i?) 
independent of luminosity and redshift (and hence out- 
side the integral). The question of whether (i?) evolves 
or not with z is still unclear but it seems likely that it 
does vary with luminosity (see Section [?]). 

We will repeat our analysis using priors centred on 
{R)^2 and 4, to show the different values inferred for 
(e). In addition we will also carry out the analysis with- 
out including obscured AGNs, using Equations [5] andfTUl 
which is equivalent to using a prior with {R)—0. 

In the case of (i?)'^2, inferred from the bright 
AGNs, we assign a pr i or ba sed on the estimate of 
iMartfnez-Sansigre et all (|2008D . This can be approxi- 
mated to a gaussian PDF with /i(_R,) =1.7 and ajf =1.2 
(their large uncertainty arising from small number statis- 
tics). 

Amongst low-luminosity AGNs, several studies suggest 
(i?)~4, but there is no explicit estimate that we can use 
as a prior. For this case, we use a ^-function prior cen- 
tred on 4. Due to the lack of scatter in (R), the posterior 
PDF for (i?)=4 obtained from the 5- function prior will 
be artificially narrower than the posterior PDF assum- 
ing fi(fi ) =1 .7 and an =1.2. This will be dealt with in 
Section 

The hard X-ray LF includes unabsorbed and Compton- 
thin absorbed AGN, so that only Compton-thick AGNs 
should be missing. The luminosity-dependent absorbed 
AGN fraction is also encoded in the LF. We do not at- 
tempt to correct for the missing Compton-thin popula- 
tion, so we are again using Equation [3l and fTOl instead of 
landini 

5.1.4. P((e)) 

For the efficiency, we reflect our ignorance by setting a 
prior in log space, P(logio[(e)]), which is flat in the range 
-3.0 <logio[(e)]< 0.0. 

5.2. The local density of SMBHs 

For p,,ioc, we assume a value of 4.61^^4x10^ Mq 
Mpc"3, following iMarconi et all ()2004l ). The hkelihood 
is then assumed to be a gaussian distribution around this 
value (Equation We compare the total p.joc to the 
energy density computed using all AGNs brighter than 
Mb < -18 (for the B-band LF) or Lx >3 x lO^^ W (hard 
X-ray LF). 

For the case of (i?)'^2, given that this value is mo- 
tivated by luminous obscured AGNs (quasars), we also 
consider only the bright end of the optical LF (AGNs 
with Mb < —23). We make the assumption that they 
are hosted by the progenitors of present-day elliptical 
galaxies, so we compare their integrated energy to the 
SMBH density in elliptical galaxies only. IMarconi et al.l 
(|2004} ) estimated the SMBH density in elliptical galaxies 
to be 70% of the total SMBH density, or 3.2+} ;^ xlO^ M© 
Mpc-3. 

5.3. The posterior PDFs for (e) 

Figure [H (panels a, b, c and d) shows the resulting 
posterior PDFs for logio[(e)], given our parametrisation. 



given our priors and given our assumed local density of 
SMBHs, p,.ioc- Panel (a) shows the results from the 
B-band LFs without any correction for obscured AGNs 
((i?)=0), using our two sets of values for C^. 

Panel (b) shows the result of using the hard X-ray LF. 
The value of u derived from the hard X-ray LF (Sec- 
tion I5.1.2P includes both X-ray unabsorbed AGNs and 
X-ray absorbed Compton-thin AGNs, but not absorbed 
Compton-thick AGNs. 

Panel (c) shows the results from the B-band LFs and 
the two different bolometric corrections, with two dif- 
ferent values of (R) (~2 and 4) used to correct for the 
obscured population (which includes the Compton-thick 
AGNs as well as the Compton-thin ones) . Note that the 
posteriors for (e) derived from {R)—0 and (i?)=4 are 
artificially narrower than that for {R)^2, due to the S- 
function priors for (R). This is dealt with in Section jSTH 

Panel d) shows the posterior PDFs for the quasars. 
It is derived from the B-band LFs, integrating only over 
AGNs brighter than Mb = —23, considering only (i?)~2, 
and comparing to the SMBH density in local elliptical 
galaxies. 

5.4. Combining different estimates for (e) 

The estimates for (e) shown in Figure 2 have been de- 
rived using different values for C^, and in the case of pan- 
nel c, different values for (R). We now wish to combine 
the estimates in each panel, to obtain our four estimates 
for (e) for each of the three scenarios: a) considering only 
unobscured quasars from the B-band LF b) considering 
the absorbed and unabsorbed AGNs from the X-ray LF, 
but not the Compton-thick ones and c) attempting to 
include all obscured and unobscured AGNs by using the 
B-band LF and correcting by {1+{R)) finally d) is the 
same as c) but only for the quasars. 

In all four scenarios, we have used two different sets 
of bolometric corrections. We have no strong reason to 
believe one set is more appropriate than the other, but 
they do yield different values of (e) . Additionally, in sce- 
nario c, have derived estimates for (e) in the presence 
of obscured AGNs based on two different priors for (R) , 
one where (i?}~2, the other with (i?)=4. Both of these 
values for {R) are believed to be accurate, but only ap- 
propriate at a certain luminosity. Instead, we have used 
a fixed value of (R) for the entire LF. We now wish to 
improve our estimate, by considering how appropriate 
each of our estimates is, and combining th em. 

We follow the arguments of lPressI (|1997| ). and consider 
our different estimates of (e) from different combinations 
of the parameters Ci, and (R) as indepedent measure- 
ments. Although these inferred values of (e) have been 
derived using the same method, and are therefore ap- 
parently not independent, the reason they differ is due 
to the different assumptions for Ci, and (i?), which are 
themselves independent. What we are doing here is to 
estimate the most likely value of (e), given the indepen- 
dent parameters assumed. 

We must now assess the appropriateness of each esti- 
mate for (e) resulting from using each combination of 
parameters. Each resulting estimate of (e), (e)^, will be 
considered a good estimate or a bad estimate, depending 
on whether the values of Ci, and (R) assumed were ap- 
propriate or not. Bad estimates of (e) are simply those 
where the standard deviation has been underestimated. 
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To each measurement (e)^, we assign the probabiUty 
that the choice of parameters was the most appropriate 
and probability (l-pi) that it was not the most appro- 
priate. We here consider all p^'s equal, simplifying to 
P- 

The posterior probability for the efficiency (e), given 
data D, can then be written as 

P(logio[(e)]|D) = 

^^^^^I^Jnp)l[^pPG. + {i-p)PBMP (12) 

where Pci and Pbi are the probability distribution func- 
tions for good and bad measurements, derived fro m ap- 
propr iate and inappropriate parameter values (see iPressI 
119971 for a derivation of this result). 

We marginalise over p. Thus, even if we do not have 
any information on the probability of a term being cor- 
rect or not, and have no objective way of assigning values 
of p, we can still proceed. We simply integrate over the 
entire range of values for p, weighted by the prior P(p). 

To estimate P(log]^Q[(e)] |D), the posterior PDF for (e) 
given our different measurements, we assume a flat prior 
for logio[(e)]. P(logiQ[(e)])/P(D) is in practice simply a 
normalisation term. For p, we assign a prior flat in the 
range < p < 1, reflecting complete ignorance. 

The good probability distribution functions Pg^ are 
those where the uncertainty is correctly estimated, so 
we use our individual posteriors for each combination of 
parameters (the dashed, dotted or dash-dotted curves in 
Figure 2 panel a, b, c or d). For the bad PDFs, Psi, we 
choose to assume they would have been correct if they 
had quoted an uncertainty twice as large. We rederive 
the curves, using exactly the same procedure and priors 
of Section O but with a flnal uncertainty which is made 
artificially to be twice as large. 

When combining the different estimates, we obtain an 
improved one that takes into account not only the differ- 
ent mean parameter values, but also the different uncer- 
tainties. For example, the potentially incorrect (under- 
estimated) width of the estimates for (e) derived from a 
(5- function prior at (i?)=4 has been dealt with. 

The resulting combined posterior PDFs for (e) can be 
seen as the solid curves in Figure [Ha, b, c and d. These 
are our best estimate for the mean radiation efficiency, 
(e), given each scenario. 

5.5. Results 

Our results for the different scenarios illustrate how 
the efficiency estimated from the Soltan argument varies 
as the population of obscured AGNs is taken into ac- 
count. In all panels of Figure 2, the solid lines show the 
result from combining the individual posterior PDFs as- 
suming different parameters. The uncertainties quoted 
below correspond to the gaussian approximation of the 
solid lines in each figure. This gaussian approximation 
is extremely accurate for the low-efficiency tail of the 
PDFs, but underestimates the high-efficiency tail typi- 
cally above the full-width half-maximum. 

When we consider only the unobscured AGNs, we find 
logio[(e}]=-1.77t°:ii (Figure 2a) or a relatively low efli- 

ciency, (e) = 1.7to;4%. This is the efficiency derived when 
no obscured AGNs are taken into account. 



Using the hard X-ray LF (Figure 2b), we find 
logio[(e)]=-1.20ig:ig or (e)=6.4i^-J%. This estimate for 
(e) includes the absorbed Compton-thin population, but 
is still missing the Compton-thick AGNs. 

Figure 2c shows the result if we integrate over 
the optical LF and consider the values of (R) for 
high-luminosity and low-luminosity AGNs, logio[(e)]=- 
l.lTio Jg or (e)=6.7tj;?%. We consider this our overaU 
best estimate for (e), because it is expected to include 
the Compton-thick AGNs and therefore to be the least 
incomplete in terms of obscured AGNs. 

Finally, when only obscured and unobscured quasars 
are considered we obtain the solid curve in Figure 2d. 
Here we are using the optical LF integrated above Mb = 
—23 and assuming (i?)~2. We are also only comparing 
the radiated energy density to the mass density of local 
SMBHs in elliptical galaxies. The results can then be 
approximated as logio[(e)]=-1.27t°:ii or {e)=bAtll%- 

The difference in efficiency between the B-band LF 
with {R)=0 and the other results shows the cosmolog- 
ical importance of obscured (or absorbed) AGNs. In 
Section [51 we discuss how this affects the estimate of 
SMBH spins, under a set of assumptions. However, from 
the different results obtained using different values of C^, 
visible in the different panels of Figure [H we can also in- 
fer that the bolometric correction is still a critical source 
of uncertainty. 

The value of (e) has been a matter of some debate, 
particularly whether its value is above or below ~10%. 
In Section 15.61 we compare to results obtained by other 
studies. We note that, from our best-estimate PDFs from 
Figure[2lD, c and d, we find P((e}<10%) ranges between 
79 and 88%. Hence, our results strongly favour values of 
(e) below 10%. 

5.6. Comparison to previous estimates 

Given the large uncertainties in (i?), in Ci, and in 
P»,ioc, our integrated approach is appropriate. Using this 
method, we obtain (e)^7%, and find ou r results consis- 
tent w ith those of iMarconi et al.l (|2004f l. Shankar ct al] 
(|2004D and lMerloni fc Heinzl (|2008D who used differential 
LFs (typically hard X-ray LFs) and local SMBH mass 
functions, and applied extra corrections to account for 
the Compton-thick population. 

Our res ults are partially in agreement with those of 
ICao fc Lil |2008), who find (e)-8% when considering the 
AGNs leading to SMBHs with masses < 10^ Mq. How- 
ever, they also find (e)^18% for > 10^ M©, whereas we 
have found that {e)^5% when comparing quasars to the 
SMBH mass density in elliptical galaxies. The SMBHs 
in ellipticals do not necessarily ma p one-to-one with th e 
SMBHs with > 10^ Mq (see e.g. iMarconi et al.ll2004 . 
but our results suggest that the mean radiation efficiency 
remains approximately constant. 

We d o not, how ever, find a good agreement with the re- 
suhs of lYu fc L u (2008), who d erive a value of (e)~16%, 
or those of IWang et al.l ()2006f ) who used a completely 
independent approach to estimate (e)=20-35%. 

Surprisingly, our res ult is also consistent with that of 
lYu fc Tremaind ()2002l ). despite the fact that they did not 
include the obscured AGN population. These authors 
used a lower value for p, Joe (2.5 ± 0.4x 10^ M© Mpc^^) 
and also used a higher value of C^, , since they did not re- 
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move the infrared component of the quasar SED. These 
factors approximately cancel the effect of including the 
term, expl aining why our estimat e of (e) is con- 
sistent with that of lYu fc Tremaind (120021 ). 

Using the intensity of the X-ray background, and as- 
suming a_t}[2i£§J^6dshift of 2 for the sources dominat- 
ing it, lElvis et all (|2002D estimated (e)> 15%. X-ray 
selected AGNs, however, have been observed to lie pre- 
dominantly at lo wer redshifts (e.g. iBarger et al.l 120031 : 
lUeda et al.ll2003f) . If a typical value of z = 0.7 - 1 is 
assumed inste ad, then the value of (e) derived following 
the method of lElvis et all (j2002f ) is reduced to -8-10%, 
in reasonable agreement with our best estimate. 

6. ESTIMATING THE MEAN SMBH SPIN 

The rotation of a black hole determines the radius of 
the innermost stable circular orbit (ISCO) and therefore 
its binding energy. The rotation is measured in terms 
of the parameter a = (in units of ^^r^), where J is 
the angular momentum, to, is the mass of the SMBH, G 
is the gravitational constant and c is the speed of light. 
The accretion efficiency (77, equal to 1 minus the binding 
energy and shown in Figure [3]) is therefore determined 
by a. 

The standard model for a thin accretion disk as- 
sumes no stresses are present between the ISCO and the 
Schwartzschild radius and neglects the effects of radia- 
tion from matter falling through this section, as well as 
assuming that any hea t dissipated in the disk i s imme- 
diately radiated away (|Novikov fc Thornill973[ ). These 
assumptions then lead to e = r\(a)^ so that the mean radi- 
ation efficiency can be used to estimate the mean SMBH 
spin. 

Relaxing these assumptions has two important conse- 
quences, e is altered and, through accretion alone the 
SMBH might reach an equilibrium spin value lower than 
a= 

When radiation inside the ISCO is considered, the 
black hole captures preferentially photons with oppo- 
site angular momentum, which limits the final spin 
to a =0.998^ (iGodfrevI [19701: iThornd [l97l . Mag- 
netic fields in the region inside t he ISCO can trans- 
port angular momentum outwards (|Blandford fc Znaiekl 
119771). possibly further l imiting the equilib rium value of 
o (jGammie et al.l l2004t iKrolik et all l2005l estimate a 
limit of -0.9^). 

The assumption that e = rj requires that no torques are 
present within the ISCO, and that the matter "plunges" 
into the black hole too quickly for any radiation to be 
significant. In this case, no energy can be extracted 
once matter has crossed the ISCO. However, the infalling 
matter can exhert a torque through magnetic fields, so 
that work is still done on the disk even after crossing 
the ISCO. The dissipation of this work can then lead to 
a higher radiativ e efficiency (e.g. iBeckwith et al.l 120081 : 
iNoble et al.l[200g l. 

In the following section, we discuss our estimates of 
(a) from (e) assuming e = 77, but from the discussion 
above these values should be considered as upper limits, 
since e is really expected to be > r]{a). For the case of 
SMBH and disk co-rotating with a - 0.9^, the ISCO 
is very close to the Schwartzschild radius, so that there is 
little time to radiate in the region inside of the ISCO. In 



addition, the gravitational redshift diminish es the energy 
of the emergent radiation. For this case, iNoble et al.l 
(|2008!) have found that e — I.O677, so the difference is 
negligible compared to our uncertainties. However, for 
lower values of a the differe nce between e and rj m ight be 
larger, even of order unity (jBeckwith et al.l[2008f l. 

Our value of (e) derived in the absence of obscured 
AGNs {{e)^1.7toi% for {R)=Q) is consistent with 
SMBHs having no rotation during periods of accretion 
((e)<5.7% so that a = 0). An important point is there- 
fore that when obscured AGNs are not included, the de- 
rived efficiency suggests that SMBHs were not spinning 
during periods of accretion. 

When we do include the obscured AGNs, we find that 
values of (e)>5.7% are likely, so that the possibility that 
the SMBHs are rotating must be considered. 

Maximally rotating black holes have a — with 
77 =42% (|Bardeenlll970[) . although in the standard ac- 
cretion disk model they are only expected to reach a = 
0.998 %^ at most, corresponding to 77 =30% (jThornd 
1 19741 ). From our posterior PDFs, P((e)>30%) ranges 
from 0.3%^ to 2%^. Thus, under our parametrisation, we 
can rule out that SMBHs are, on average, maximally ro- 
tating with >98% confidence. More recent work suggests 
the upper limit on a might be lower, —0.9 ^^^, and with 
a radiative efficiency ~15% (Noble et al.''2008'). Again, 
from our posterior PDFs, P((e}>15%) ranges from 2% to 
8%, so we rule out (a) > 0.9^?^ at the >92% confidence. 

For spinning black holes, (77) will vary depending on 
whether the ISCO is co-rotating or counter-rotating with 
respect to the black hole spin (see Figure [3]). To estimate 
the mean value of a two possible limiting scenarios can 
be envisaged: 

A ) The SMBHs grow from the seed in one long episode 
of accretion. In this case one single accretion disk is ex- 
pected to be able to provide enough mass to turn the 
seed into a SMBH. The angular momentum of the seed 
is negligible compared to that of the accretion disk, so 
the disk will inevitably force the black hole to co-align 
(jScheucr & Fciler 1996; King et al. 2005). It is therefore 
safe to estimate the efficiency assuming the black hole 
and ISCO are always co-rotating. In this case, our best 
estimate of {e)—6.7^\'^% corresponds to SMBHs rotat- 
ing with (a) — 0.25to'25^^T^ during periods of accretion 
(Figure [3l solid line). 

B) The SMBHs grow in many shorter periods of ac- 
cretion. The implication here is that the angular mo- 
menta of the accretion disk and black hole are compara- 
ble and are also initially randomly oriented with respect 
to each other. Depending on the initial configuration 
and ratio of angular momenta, the bl ack hole and disk 
will either co-align or counter-ahgn (jKing et al.l 120051 : 
iLodato fc Pringld I2OO6I ). It is expected that co- and 
counter-aligned accretion is equally probable, so when es- 
timating a, a mean {ri{a)) from co- and counter-rotating 
or bits should be use d (the symmetrised curve suggested 
by iKing et all l2008l shown in Figure [3] by the dashed 
line). In this case, for our best estimate of (e) we inferr 

(a) = 0.59^^59^- 

These probabilities have been calculated from the actual PDFs, 
not the gaussian approximation which underestimates the proba- 
bility in the tails of the PDFs. 



8 



A. Martmez-Sansigre & A. M. Taylor 



Mergers involving SMBHs are expected to prevent 
a from reaching the high est values (see the work of 
iHughes fc Blandfordl[2003l . who find half of their param- 
eter space will lead to a < 0.4 '^™* ). For both accretion 
scenarios, and for a ^ 0.8^^5^, the dependence of r/ on 
a is weak so that it is difficult to constrain (a) from our 
estimate of (e). However, we have found that maximally 
rotating SMBHs are not favoured since P((e)>30%)<2%, 
and our best estim ates for a are in reasonable agreement 
with the results of iHughes fc BlandfordI (|2003[ ). 

Our estimates for (a) are nevertheless also consistent 
with the value of a > 0.5^^ es timated for the SM BH in 
the Milky Way, Sagittarius A* (jGenzel et al.ll2003l who 
explained the short near- infrared flare period of 17 min- 
utes as being close the ISCO for co-rotating disk and 
SMBH). A SMBH with a > 0.5^ and a prograde or- 
bit has an accretion efficiency rj >8%. Thus, the limit 
on a inferred for the individual SMBH Sagittarius A* is 
in good agreement with our best estimate of the mean 
radiation efficiency, (e), assuming the standard thin disk 
model and ignoring the magnetic fields. 

Indeed, we remind the reader that we have estimated 
(a) assuming e — -q, but at the beginning of this Sec- 
tion we have argued for e > rj. Hence, our estimates of 
(a) can be consider ed upper limits, since for a given true 
value of g, e > 17(a) ([Krolik et al.ll2005HNoble et alll2008l : 
iBeckwith eraIll2008D . For a ~ 0.9^t*-, it seems that e is 
only slightly larger than 77, but for lower values of a the 
difference is expected to be larger, since the time avail- 
able for radiation is longer and the gravitational redshift 
is less severe. This only strengthens our main conclusion 
that our inferred value of (e) does not require accreting 
SMBHs to be rotating particularly rapidly, on average, 
during periods of accretion. 

7. THE MEAN E-FOLDING TIME 

Given our assumption that AGNs grow by accretion, 
their mass will grow exponentially provided that the lu- 
minosity is proportional to the mass. Thus, if AGNs ra- 
diate at a fraction A of their Eddington luminosity, then 
they will grow exponentially with an e-folding time: 



A(l - e) 47rG'mp ^ ' 

where rrip is the mass of the proton and gt is the Thom- 
son electron scattering cross-section. Our estimate of (e) 
is completely independent of A, and it is only to estimate 
T that we require an estimate of the Eddington ratio A. 

Estimates of A exist for unobscured quasars, derived 
from their bolometric lumino sities and virial-ma s s esti - 
mates for the SMBH mass. iMcLure fc Dunlod (pOOl ) 
found A 15-0.5 in the range 0.1 < z < 2.1, 
with significant scatter. We approximate the distribu- 
tion of logio[A] to a gaussian with (Ltiogjj,[A]=-0.6 and 

0'logio(A)=0.4. ^ 

When combined with our estimate of (e) for quasars 
(5.4l^;2%), this leads to an estimate of logio[(T)/Myr] = 
2.0 ± 0.4, or (t) = lOOtgo^ Myr. This value is consis- 
tent with (although at the high end of) other current 
constraints from quasar clustering measurements, from 
arguments based on coevolution of quasars, SMBHs and 
host galaxies, from arguments of simultaneous trigger- 



ing of quasars and starbursts and from the lengths of 
radio jets. Thes e estimates al l suggest r is in the range 
10-100 Myr fsee lMartinill2004 for a review). 

Such a long value of t suggests SMBHs require a long 
time to build up their mass. At z = 6.4, SMBH masses 
of 2-6x10^ MfTi h ave been estimated (Willott et al..,,2003l : 
iBarth et al.ll2003l ) , when the Universe was only 840 Myrs 
old. This suggests that the seeds that formed the z = 6 
quasars might have been as large as ^ 10^ — 10^ M©. 

The estimates of the SMBH mass and particularly our 
estimate of the e-folding time have very large uncertain- 
ties and the quasars at z — 6.4 could have shorter e- 
folding times. They might be accreting closer to the Ed- 
dington li mit than A^O.25, the me an of the 0.1 < z < 2.1 
quasars (jMcLure fc DunlaDll2004D . or radiating less effi- 
ciently than ~7% . This would then imply significantly 
smaller seeds. For example, if r ^60 Myr (40 Myr), the 
seeds could be -10^ M© (~1 Mq). 

8. SUMMARY 

We have revisited the mean radiation efficiency from 
accretion onto SMBHs, (e), quantifying the effect intro- 
duced by taking into account a large population of ob- 
scured AGNs. We estimate the cumulative energy den- 
sity of the AGN population from the unobscured AGN 
LF by applying appropriate bolometric corrections and 
correcting for the obscured population. The efficiency 
(e) is estimated by comparing this energy density to the 
SMBH mass density in the local Universe. 

When no obscured AGNs are included, we find 
{€) = 1.7toi% which does not require SMBHs to rotate 
at all during periods of accretion ((«)= 0). 

Using the hard X-ray LF, which includes absorbed 
Compton-thin AGNs, we find (e)=6.4t^j%. Using 
the B-band LFs, correcting for the obscured popula- 
tion (which includes the Compton-thick AGNs), we find 
(e)=6.7;ti;i%, which we consider to be our best estimate. 

We have also derived (e) for the powerful AGNs only 
(the quasars), and find {e}=5.4:'tl'^%- 

Depending on whether the SMBHs and accretion disks 
are assumed to be always co-rotating, or instead to 
counter-rotate in half of the cases, we estimate the ro- 
tation parameter during periods of accretion to be (a) = 
0.25t|];30^^ or (a) = Q.59t°oil^, respectively 

Efficiencies >30% (>15%) are ruled out at >98% 
(>92%) confidence, so the SMBHs are not, on average, 
maximally rotating during the periods of accretion. Effi- 
ciencies <10% are favoured, with 79%-88% of the prob- 
ability density in this region. 

Finally, combining our best estimate for (e) with es- 
timates of the Eddington ratio for unobscured quasars 
yields an estimate of the e-folding time of (r) = 
lOOtgo^ This suggests that the seeds of z = 6 

quasars might have been very large, potentially as large 
as ~ 10^ Mq. However, the large uncertainties in r trans- 
late into order-of-magnitude uncertainties for the masses 
of the seeds, which could be as small as ~ 1 M©. 

We are grateful to John Miller, Andrew King and 
Giuseppe Lodato for useful discussions, to Ross McLure 
for access to the distribution of Eddington ratios, to 
Fabio Fontanot and Mark Sargent for comments on the 
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Fig. 1. — Cumulative energy densities for the optical LFs (left) and hard X-ray L F (right). We r emove any dependence on Cu by dividing 
them by the maximum value (at z = 0). The maximum values are taken from the ICroom et all II2004I) LF (left) and from the lUeda et al.l 
I I2OO 3) LF (right). The Wolf et al] II2003I ) LF, defined only for ^ > 1, agrees well with the other two LFs assuming pure density evolution 
(PDE). However, pure luminosity evolution (PLE) predicts a somewhat larger value for Ub- The optical LFs include only unobscured 
AGNs, while the X-ray LF also includes Compton-thin absorbed AGNs. 



Unobscured AGKs only, 

(a) 




Inc. compt.on — thin AGNs, 11^^ 

(b) 




Fig. 2. — Posterior PDF for (e) derived from the optical LFs for unobscured AGNs only (a) and using the hard X-ray LF, wit h no 
correction for Compton-thick absor bed AGNs (b). The posteriors derived from the bolometric corrections of IMarconi et aLl 120041 ) are 
marked as 'M04', whereas those from lHopkins "etaLl 1 I2OO7I) are marked as 'H07'. The combined posteriors, derived in Section l5.4l are shown 
as the solid curve. For panel (a), the optical LF is integrated over all mag nitudcs brighter than Mb = -18 {L^b>Sx10'^° W Hz"!), 
whereas for panel (b) the LF has been integrated over all Lx > 3 X 10^* W. 



Inc. obscured AGNs, (l + (R)}f/p Inc. obscured quasars, (l + (R))f/, 




Fig. 2. — Continued. Posterior PDFs for (e) derived from the optical LFs correcting for the obscured AGNs using (R). In panel (c), 
the optical LF is integrated over all magnitudes brighter than Mb = —18, and two different values of {R) are used ((i?)~2 and =4). We 
consider this our best estimate for (e). In panel (d), the LF has been integrated over all magnitudes brighter than Mb = —23, and only 
(R)~2 is used, to estimate the value of (e) from the bright AGN population (the quasars). 
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Fig. 3. — Accretion efficiency ?7 as a function of rotation parameter a. The solid curve represents the case of the SMBH co-rotating with 
the inner part of the accretion disk, so the innermost stable circular stable orbit (ISCO) is prograde. The dotted line represents the case 
of counter-rotating SMBH and inner accretion disk (the ISCO is retrograde). The dashed curve represents the mean of both curves (the 
"symmetrised" curve, see IKine et al.|[2008) . 



